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We study the density of states and magnetotransport properties of disordered Weyl semimetals, 
focusing on the case of a strong long-range disorder. To calculate the disorder-averaged density of 
states close to nodal points, we treat exactly the long-range random potential fluctuations produced 
by charged impurities, while the short-range component of disorder potential is included systemat¬ 
ically and controllably with the help of a diagram technique. We find that for energies close to the 
degeneracy point, long-range potential fluctuations lead to a finite density of states. In the context 
of transport, we discuss that a self-consistent theory of screening in magnetic field may conceivably 
lead to non-monotonic low-field magnetoresistance. 


I. INTRODUCTION 

Coulomb disorder plays dominant role in determin¬ 
ing thermodynamic and kinetic properties of doped 
semiconductors due to its long-range nature.® This 
physics becomes especially intriguing in the context of 
Dirac materials, namely solid state systems that host 
quasiparticles with linear dispersion near the band¬ 
touching degeneracy points, such as that in (i-wave 
superconductors,® gr^hene,®^ topological insulator^^ 
and Weyl semimetals.®® While the former three classes 
of systems have been intensively studied over almost a 
decade by now, materials that can be identified as Weyl 
semimetals have been discovered only fairly recent ly.^^^HIll 

Theoretical studies of the effect of disorder on a single 
Weyl node have a long history,^^ and thi^roblem has 
been revisited in a number of recent works.^^^^^The pre¬ 
vailing point of view is that a weak disorder has negligible 
effect on the density of states, which vanishes quadrati- 
caly with the energy counted from the nodal point. This 
behavior persists up to a certain critical disorder strength 
beyond which the density of states acquires a finite value 
at zero energy. Obviously, such scenario should have se¬ 
rious implications for the transport properties of WSMs 
that rely on their semimetallic nature (a vanishing den¬ 
sity of states at the nodal points).I^^M^ At strong disorder 
singular transport features are eliminated. It is thus im¬ 
portant to have a controlled theory of disorder effects 
on a single nodal point. Conventional treatments based 
on the self-consistent Born approximation, though can 
be qualitati yely correct, are uncontrolled near Dirac or 
Weyl points.l^ 

Magnetotransport properties of WSMs are also ex¬ 
pected to have a highly unusual character owing to 
the unconventional Landau quantization for the three- 
dimensional Dirac spectrum, and the multiband na¬ 
ture of Weyl semimetals. In particular, a strong 
classical negative magnetoresistance was predicte d for 
Weyl systems with long intravalley relaxation times,^^^^^ 
as well as robust linear positive magnetoresistance in 
the quantum limit, when a single Landau level is 


occupied.l^ Both ty pes o f behavior were recently ob¬ 
served in experiment on a candidate Dirac (“double- 
Weyl”) material Cd 3 As 2 , yet at this point it cannot be 
stated with confidence whether the experimental obser¬ 
vations are in full correspondence with theoretical mod¬ 
els. For various recent theoretical and experimental as¬ 
pects of magnetotransport results in Dirac semimetals 
see Refs. I37ll46l 

In light of the recent experimental advances that trig¬ 
gered a flood of theoretical work, our essential motivation 
is to investigate the role of Coulomb impurities on prop¬ 
erties of WSMs. In contrast to the model of short-range 
disorder, which so far received most of theoretical at¬ 
tention, we obtain results that are controlled by a small 
parameter in the limit of strong disorder. We also ob¬ 
serve that due to the delicate interplay of Landau quan¬ 
tization and screening effects, the low-field magnetore¬ 
sistance may display nonmonotonic behavior. The latter 
feature is highly sensitive to the strength of interaction 
and thus is not universal. 


II. MODEL OF A WEYL SEMIMETAL WITH 
LONG-RANGE DISORDER 

The quasiparticle dispersion in the vicinity of a single 
isotropic nodal point can be described by the effective 
Hamiltonian { h = c = 1 throughout) 

H = va-pPu{r), (I) 

where v is the Fermi velocity, and a and p are vector 
of the Pauli matrices, and momentum measured rela¬ 
tive to the nodal point. The total potential as seen by 
an electron at position r due to randomly distributed 
charged impurities at positions Vi is taken in the form of 
a Yukawa-like potential 

E g2 

0(r —rA, (fir) = —exp(—/^r). (2) 

i 

Here index i labels impurities, e is the dielectric constant 
of the host material, and n is the inverse screening ra¬ 
dius. We assume that impurity positions are completely 
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random, and calculate the correlator of the disorder po¬ 
tential, 


{u{r)u{r')) 


/47re^\^ f d^q nex.p{iq{r — r')) 

\ € ) J (27r)3 (g2 + ^2)2 


27rne^ . , 

— 5 — exp(-/^ r 


r'\). 


( 3 ) 


where n is the concentration of impurities. As follows 
from Eq. © , the correlation length of the disorder poten¬ 
tial is given by the screening radius, and the magnitude 
of its fluctuations is given by {v?{r)) = = 2'Kne^ j k. 

The Gaussian approximation is sufficient to describe the 
disorder potential fluctuations when there are many im¬ 
purities in the screening volume, ^ 1. The model 

of a disordered Weyl semimetal described above will be 
used in the forthcoming calculations. The validity of the 
Gaussian approximation will be checked a posteriori^ by 
considering the screening of Goulomb impurities by elec¬ 
trons moving in a self-consistently screened disorder po¬ 
tential, see Eq. (33), and the text below it. A detailed 
account of self-consistent screening in the present prob¬ 
lem is given in Ref. [23l 


III. DENSITY OF STATES AT THE NODAL 
POINT 

In what follows we consider the density of states (DoS) 
close to the nodal point, 5^0, where 5 is the electron 
energy. 

DoS at the nodal point crucially depends on the re¬ 
lation between the strength of disorder, and the energy 
scale associated with its finite correlation length, kv. For 
w hzv, the disorder is weak, and the Born approx¬ 
imation can be used, with the expected result of zero 
DoS at the nodal point. For strong disorder, w ^ kv, 
self-consistent Born approximation predicts a finite DoS 
at the noda l p ointbut the approximation itself is 
uncontrolled 

To obtain systematic results, controlled by a small pa¬ 
rameter, for the DoS and later for transport properties, 
we restructure the perturbation theory for the disorder 
potential similarly to the way it was done for heavily 
doped semiconduct or s.flEll A crucial observation is that 
for w ^ Kv, the landscape of the disorder potential is 
smooth, in the sense that the typical length scale at which 
it changes, 1//^, is large compared to the typical electron 
wavelength A ^ v/w. (The latter estimate assumes that 
it is the total energy of the electron that is close to the 
nodal point, thus the kinetic energy is of the order of 
w.) This implies that the electron motion is semiclassi- 
cal. Therefore, one may use the Thomas-Fermi approxi¬ 
mation (rather than the clean system) as the “reference 
point” around which a perturbative expansion is devel¬ 
oped. 

To illustrate the use of the Thomas-Fermi (TF) ap¬ 
proximation in this case, we consider the single particle 


Green’s function for the Hamiltonian 0 : 

[n(r)] =-^, (4) 

e± — va p — u[r) 

where e± = 5 ± iO, and u{r) plays a role of the effective 
local chemical potential. In general, since the electron 
momentum is an operator, one cannot simply average the 
Green’s function over the distribution of u{r). However, 
in the semiclassical limit we can treat u{r) as a number 
within TF-approximation and thus obtain 

= [ duF[u] - - -, (5) 

J e± — va ■ p — u 

where the distribution function for the local values of the 
disorder potential, F[u], is defined according to 

F[u\ = {6{u — u{r))). (6) 

Taking the imaginary part of the retarded component of 
(G) gives the density of states, 

r'is) = J duF[u]uQ{s — u), ( 7 ) 

where r'o(^) is the DoS of the clean system. In the case 
of a Gaussian disorder potential defined in Eq. (|^, one 
obtains 


^ exp{-u^/2w^). 

yZiTW 


(8) 


Using i^o{e) = g f - vp) = ge‘^/2TT‘^v^, with g 

being the degeneracy of the nodal point for a given Weyl 
system, and performing an elementary integration, one 
finds 


= (9) 

that is, a nonzero value for e = 0. 

In the context of doped semiconductors, the semiclas¬ 
sical equation for the density of states, Eq. 0 , was sug¬ 
gested long ago, e.g. in Refs. 113® and has been re¬ 
cently used to describe DoS and nonlinear screening in 
Dirac systems in Ref. 

We motivated equation 0 on physical grounds, but it 
can be rigorously shown that it is exact in the vn/w 0 
limit. One can formally send k ^ 0 while keeping w 
constant by considering ^ 0, n ^ oo limit with 
kept constant (see below). In this limit the disorder cor¬ 
relator is {u{r)u{r')) ~ and describes disorder with 
infinite correlation length. The problem of finding the 
single-particle Green’s function for such a disorder poten¬ 
tial (the so-called “Keldysh model”) is exactly solvable, 
since the impurity lines in the standard technique do not 
transfer momentum, thus all diagrams to a given order 
have the same value. Fig. The full resummation of 
the perturbation theory becomes a combinatorial prob¬ 
lem, solving whiclPIl requires determining the number of 
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{ w‘^, real space 

(27r)^5(q)w;^, momentum space 


FIG. 1: Diagrams that correspond to the second-order disor¬ 
der correction to the Green’s function. For the impurity line 
that does not depend on spatial coordinates, {u{r)u{r')) = 
all the diagrams have the same value of w^Gq, where Go 
is the unperturbed Green’s function. 


ways n impurity lines can be attached to 2n H- 1 bare 
Green’s functions in the order of perturbation theory 
(the second order diagrams are shown in Fig. [^. The 
results for the disorder-averaged single-particle Green’s 
functions and the DoS are precisely the above Eqs. (§ 
and 0. 

In what follows we consider the corrections to the result 
of Eq. following the treatment of Ref. 1471 


A. Improved perturbation theory for long-range 
disorder 


To build an improved perturbation theory for the po¬ 
tential defined by Eq. ([^, we formally rewrite the latter 
equation as 

{u{r)u{r')) = w‘^ (exp(—/^|r — r'\) — 1). (10) 

If the second term in the right hand side of the above 
equation is neglected (the corresponding condition to be 
formulated below), one recovers the case of an infinite- 
range disorder potential, which was discussed above. 

Going beyond the Keldysh model, when is it appro¬ 
priate to treat the second term in the right hand side 
of Eq. (10) as a small correction? It is clear that one 


should require that K\r — r'\ 1 for the relevant values 

of |r — r'|. Since the solution for the semiclassical prob¬ 
lem is represented by an ensemble average over samples 
having a random, but spatially uniform chemical poten¬ 
tial of order re, see Eqs. (§ and ([^, the relevant value of 
\r — r'\ is of the order of the wavelength corresponding to 
the energy re, \r—r'\ ^ vfw. Hence the corrections to the 
semiclassical approximations are expected to be small for 
Kv/w 1, coinciding with the condition for the validity 
of the semiclassical approach formulated above. 

Our next immediate goal is to find the leading ^ 
0{vn/w) correction to Eq. © due to the short-range 
component of disorder potential. As we are going to see, 
it is regular and free from divergencies near the Dirac 
point. The full disorder-averaged Green’s function can 
be written 


(^qR(A)) ^ J 


1 


e± — u — va ■ p — H 


(11) 


Here, for every value of the self-energy E is determined 
from the usual diagram technique for disordered metals, 
in which the role of the bare Green’s function is played 
by Eq. 0 with u{r) replaced with and the impurity 
line (in momentum space), D(gr), is determined only by 
the second term in the right hand side of Eq. ( [Tq| ) : 

J d^qD{q) = 0. 

In Eq. 5{x) is the Dirac delta-function. Its appear¬ 
ance signals that the long-range part of the impurity po¬ 
tential is subtracted from D{q), as it has been already 
treated exactly by the averaging over the values of u. 
We will see below that after such removal, a simple Born 
approximation is valid for the short-range part of the dis¬ 
order potential if kv w. 

Nonetheless, below we set up the self-consistent Born 
approximation (SCBA) for the self-energy for the modi¬ 
fied disorder correlator, Eq. ( p!^ . We will not go beyond 
the first Born approximation in actual calculations, but 
having the full SCBA set up will allow us to contrast 
the present treatment with the conventional SCBA for 
long-range disorder. 

Using the symmetry-suggested decomposition for the 
self-energy 


E — Eq T E]^tT • TXp^ 


(13) 


where rip is the unit vector in the direction of p, and 
introducing = s — u±i0^ for the retarded component 
we have 


2 poo 

^o{p) = - dqq^ / dx 

TT Jo J _i 


Kw'^go{q) 


-w^goip), 

2 


poo p-\- 

^i{p) = - dqq^ / 
Jo J -1 

-w'^giip), 

where we have introduced 


OO p-\-l 

^ ' dx 


(p^ — 2pqxY 

(14) 

i<iw‘^xgi{q) 

{p^ P q^ P n? — 2pqxy 

(15) 


go 


9i 


So 


(e+-Eo)2-(^pH-Ei)2’ 

vpP El 

(e+-Eo)2-(i;p + Ei)2- 


(16) 

(17) 


To simplify equations, we first introduce new functions 
Xp = — Eq and Yp = vp-\-Yi. Second, we complete x- 

integration in Eqs. (14) and (15) and observe that upon 
changing q ^ —q the corresponding integrands are even 
and odd functions respectively. We thus extend the p- 
integration to the whole axis such that Xp is considered 
as an even function while Yp is odd. As a result we have 
to solve the following set of coupled integral equations 




“7-00 


KW 

Yp = vp^ - 

TT 


— CX3 

-|-oo 


Ko{p,q)Xgdq 

X^-Y^ 

q q 

Ki{p,q)Yqdq _^2 
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where the kernels are given explicitly by 


Ko{p, q) 

Ki{p,q) 


q 1 

p {p — q)^ + ’ 


1 1 , 
+ -— In 
Apq 


ip - qY + 


/ jp - qf + K^ \ 

\ip + q)'^ + K? ) 

(19) 


At ^ 0 limit within the self-consistent Born 
approximation and Keldysh model 


It is instructive to understand how the limit of very 
long correlation length, ^ 0,re = const, is reached 
within various approximation schemes. Firstly, the re¬ 
sults for the Keldysh model, Eqs. (§ and 0 , are exact 
in this limit, as we have already discussed. In particular, 
this means that the quantities X and K, that determine 
the self-energy at the SCBA level for the modified disor¬ 
der correlator, must attain bare values Xp = = vp 

when Eqs. (18) are solved. That this is indeed the case 


can be seen noting that in the /i: ^ 0 limit the kernels 
Ko,i satisfy 


lim =^iP- 

K^O TT 


q)- 


( 20 ) 


This implies that the second terms on the right hand sides 
of Eqs. (18) are exactly cancelled by the third terms. In 


fact, all self-energy corrections for the modified disorder 
correlator, not just the SCBA ones, vanish in the k ^ 0 
limit: This is the essence of the Keldysh model, in which 
the full diagrammatic series is summed up exactly i n the 
^ 0 limit, to yield Eq. ^ for the Green’s function.^^^Hl 
We can contrast the Keldysh model with the con¬ 
ventional SCBA, in which the full disorder correlator, 
Eq. is used as the impurity line. This approxima¬ 
tion can easily be obtained from Eqs. (18) by setting 


^ = e — lA ^ e, and dropping the third terms on the right 
hand sides of the equations for X and Y. To distinguish 
this case from the rest of the paper, we denote X and Y 
for SCBA with full disorder correlator as X^ and Y^. 
To obtain simple analytic results, we limit ourselves to 
e = 0. In the k 0 limit the SCBA equations are 


= i5- 


Fp =vp-\-ur 




(XB)2-(F/)2’ 

Yp 

[x^Y-iY^y 


( 21 ) 


The infinitesimal S is kept as a reminder that ^Xp > 0. 
Xp and Y^ can be easily found noting that the former 
is purely imaginary, while the latter is real: 


Xp = i0{2w — vp)y w‘^ — 


{vpY 


= y + ^{vp - 2u;) 


(vpy 


— 


( 22 ) 
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FIG. 2: Full spectral weight at zero energy in the k ^ 0, w = 
const limit: Solid line - Keldysh model, dashed line - SCBA 
with full disorder correlator. 


The corresponding DoS can be found as 

(Yp 


i/-(e = 0) = -^S J 




(27r)3 (XS)2-(yB)2’ 


which upon using Eqs. (22) gives 


^^{e = 0) = 


w 




(23) 


(24) 


In this treatment inter-valley scattering is neglected, and 
the result Eq. (24) should be understood as the DoS per 
node. 

We see that the SCBA with the full disorder correlator 
reproduces the exact answer for DoS in k ^ 0 limit. This 
means that if one naively tries to correct SCBA result for 
DoS by including the low-order interference corrections 
to the self-energy,!^ the result will become worse, not 
better. 

What the conventional SCBA does not capture is the 
form of the spectral weight that results from summa¬ 
tion of the full diagrammatic series in the k 0 limit. 
As a function of momentum for a given energy e, the 
spectral weight for the Green’s function of Eq. 0 is 
a Gaussian peak centered at \e/v\^ having a width of 
w/v. This exact result is clearly not reproduced by the 
conventional SCBA at low energies, Eq. (22). Eor il¬ 


lustration purposes, we plot the total spectral weight 
(for conduction and valence bands) at zero energy en¬ 
ergy for the retarded Green’s function within the Keldysh 
model, Eq. and for the Green’s function in the ( 


tional SCBA. The latter is obtained using Eqs. (22) as 

G^(e = 0,p) = - Y^arip)-^. 


C. Leading correction to the semiclassical density 
of states 


In this Section we calculate the leading order correc¬ 
tion to the semiclassical DoS, Eq. 0- To this end we 
use Eqs. (18) to find the first perturbative correction to 
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Xp and Yp^ that is, restric ourselves to the first Born ap¬ 
proximation. It is sufficient to find the leading 0{vk/w) 
term in the DoS. To zero order in w we take Xp = and 
Yp = vp and seek the first order correction by iterations: 
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- {vqy 
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vp 




- {vpf 


,(25) 

(26) 


After performing the integrations and after some further 
algebra we find 


5Xp = w 


i+ 


i+ 


5Yj,= 


— {vpY (C+ + ivnY — {vpY 
w‘^ r {vpY + {vtiY — {vpY 


vp 


iw^vn 


2(vpY 


+ ivnY - 

+ ivK, ■ 


In 


{vpY 

- vp^ 


a 


{vpf 


+ ivn — vp 


(27) 


(28) 


These expressions appear very singular, but one must 
bear in mind that it is only the full disorder-averaged 
Green’s function that determines the DoS, and all singu¬ 
larities of the present perturbation theory are removed by 
Gaussian integration over u in the expression for the re¬ 
tarded Green’s function, Eq. (11), expanded in the usual 
series in powers of E. 

Knowledge of 5Xp and 5Yp allows one to find correc¬ 
tion to the density of states. We are only interested in 
the leading Su ^ 0{kv/w) correction. Explicit calcula¬ 
tion shows that the only term contributing to the DoS 
correction is 6iy oc ^"^pPogi^^p- Eollowing the standard 
path (expanding the Green’s function to linear order in 
K, taking its imaginary part, and performing the final 
momentum integrations), one finds 


Si'{s) = 


w 


VK, 


Att^v^ s — u iO 


(29) 


Einally, averaging this result with F[u] from Eq. ^ we 
arrive at 

(30) 

Thus the first perturbative correction to the DoS is regu¬ 
lar and controlled by the smallness of VKjw^ as expected 
from the preceding discussion. We note that at no point 
did we have to introduce an artificial momentum cut-off, 
as is commonly done while treating the system in the self- 
consistent Born approximation. Taking into account the 
second Born approximation, or the first interference cor¬ 
rection to the self-energy all lead to DoS-corrections that 
involve at least one additional power of VKjw compared 
to Eq. (30), and can be ignored in the current treatment. 


In order to establish the limits of validity of our results 
we need to determine the typical value of w. This can 
be done with the help of self-consistent Thomas-Eermi 


approximation. Eirst recall that k = J^Ke^vjt and 
vj^ = 2TTne ^!which follows from Eq. Using now 
Eq. for the density of states we find an equation that 
defines the typical magnitude of the disorder potential 
fluctuation as a function of density n in the form 

, (31) 

9 


where a = /ve is the effective interaction parameter. 
The most interesting case for us is the vicinity of the 
“neutrality” point where 5 <C re, so that 




(32) 


which sets the typical scale of the density of states and 
inverse screening length, respectively, 

uc^ain?^^lv, (33) 

As discussed above, the Gaussian model of disorder 
amounts to a condition ^ 1/a^ ^ 1. Pertur¬ 

bative treatment of the short-range disorder requires 
vn/w ^ yW ^ 1. Both of these conditions can be 
satisfied provided that a < 1. Einally, we discuss va¬ 
lidity of the semiclassical approximation that requires 
\d\/dx\ ^ 1. The latter can be estimated as follows, 
dX/dx ^ {v/v?){du/dx) ^ {v/w‘^){w/ ^ yW ^ 1. 


D. DoS in weak magnetic fields 


This semiclassical formalism can be readily extended 
to the case of finite magnetic field. The generic analysis 
of the density of states is quite involved. Here we will be 
primarily interested in the regime of weak fields, when 
magnetic length Ib = is large compared to the 

correlation radius of the disorder potential Ib^ ^ In 
this case we can use the same approach we have used 
for the magnetic-field-free case. In the strong-field limit, 
Ib ^ the treatment of Ref. |52] should be used. (See 
also Ref. |53 for an alternative quasiclassical approach to 
transport in the strong-field limit.) 

The spectrum of Landau levels for the three- 
dimensional Weyl equation has the form 


{Pz ) 


£o{Pz) 


± 

VPz^ 


+ (vp^P, 


TO = 1,2... (34) 
(35) 


where we chose a particular chirality for the zeroth Lan¬ 
dau level without loss of generality. We concentrate on 
the lowest energy limit where 


^(0) = 


9 

27r/^ 7-00 


dpz 

27r 


+ 00 

^ ^ d{£m{Pz) ~ '^)- 

m= —oo 


(36) 


Averaging this expression over u with the help of Eq. (§, 
followed by the Gaussian integration over momentum pz 
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and the geometrical series summation over m one finds 

^(0) = ^ - coth (. (37) 

{2tiIbYv \2l%w^) ^ ^ 

Expanding this result in the low field limit one recov¬ 
ers Eq. ® to the lowest order with a positive correction 
term Su/u{0) = (v/Ib)^/I 2w^ oc This finite field 
correction ultimately renders the change in the screen¬ 
ing radius of the disorder potential, which in the case 
of Coulomb impurities will have an important effect on 
magnetoconductivity. 


IV. TRANSPORT 


In this section we briefly discuss what self-consistent 
theory of screening in WSMs gives for the low- 
temperature transport, namely we concentrate on the 
regime of T ^ In the semiclassical limit, pro¬ 

vided that the mean free path I for electron scattering is 
large, kpl ^ 1, where kp is the Eermi momentum, the 
conductivity can be calculated from the Drude formula 
with energy dependent scattering time r(5), density of 
states z^(5), and cyclotron frequency 


f ds v‘^u{s)r{s) 

67r J AT cosh^ (^) 1 + (^) * 


(38) 


The ensemble-averaged transport scattering time for 
Coulomb impurities within the Boltzmann approxima¬ 
tion can be calculated using the full disorder correlator, 
Eq. since the 6{q) part of D{q) from Eq. cannot 
lead to momentum transfer. The transport scattering 
time is thus given by 


r ^ = 


Tiun 


k% 


[ (i6>|0(2/ci? sin^/2)|^ sin6>(l — cos^ ^) 

Jo 

[(2a2 + l)ln(l-ha-2)-2], (39) 


where a = K/2kB and 0(g) is the Eourier transform of 
Eq. 0 taken at momentum q = 2kF sin^/2, whereas the 
semiclassical expression for the cyclotron frequency for 
the linear spectrum is 


ujc = v'^/{l%e). 


(40) 


According to our earlier estimates of the screening radius 
one immediately concludes that n/kp ^ \/^ ^ so it 
is sufficient to retain only the large logarithmic factor in 
Eq. (39). Then, the zero field conductivity is estimated 


(up to an overall numerical factor) as follows 


(j e^V^u{u)T{u) 


2 ^ 1/3 


e n 


ln(l/a)’ 


(41) 


where we took the typical value of potential u ^ w and 
used Eq. (32). One can check that derived expression 


indeed corresponds to the Boltzmann limit. We note in 
passing that a local conductivity tensor is appropriate 
for the situation considered here, since Eqs. (32), (33), 


and (39) imply that vtk. ^ 1/| Inaj ^ 1; in other words, 
the transport mean free path is much shorter than the 
screening length. Eurther, it is worth mentioning that 
the dependence of a on the interaction parameter a was 
a subject of controversy since different results exist in the 
literature. 

The field dependence of the conductivity comes from 
two factors. The first is obviously due to the cyclotron 
motion, which we estimate by expanding the denomina¬ 
tor of Eq. 




~ -[Wc{u)T{u)f ~ - 


1 


1 




(42) 


Curiously, there is another contribution, of the oppo¬ 
site sign, that stems from the already mentioned fact, 
that magnetic field modifies the screening radius. Since 
the product vr is inversely proportional to hi{kF/ti) and 
K (X y/i) then it follows from the low field asymptotic of 
Eq. ([37[) that for this mechanism 


Sa{B) Sk{u) 


1 


1 


ln(l/ a) a‘^ 1^(1/ 


(43) 


Eor weak interactions, a ^ 1, Eq. ( |42| ) always dom¬ 
inates by a parametrically large factor and leads to 
positive magnetoresistance at weak field that eventu¬ 
ally crosses over to qu antum regime at higher fields, 
o-{B) ^ However, for some material sys¬ 

tems, where interaction is moderately strong, a ^ 1, one 
in principle could envision a scenario where the fate of 
the sign of magnetoconductivity will be determined by 
the overall numerical prefactors in Eqs. (42) and (43). 


Obviously the criterion ol ^ 1 limits the applicability of 
our estimates. 

An elaborate theory of magnetoresistance in WSMs 
was recently developed in Ref. 0^ in various regimes de¬ 
termined by the relation between the relevant energy 
scales T, v/Ib^ and Our estimates for Sa{B) at 

weak field are consistent with that of Ref. |46] with an 
addition of the contribution governed by Eq. (4^, where 
in Eqs. (41)-(43) we also retained functional dependence 


on the interaction parameter a. 


V. SUMMARY 

To conclude, we have studied the density of states and 
transport properties of Weyls semimetals with the long- 
range Coulomb disorder. Provided that the semiclassical 
limit is satisfied, DoS induced by the long-range disor¬ 
der potential fluctuations can be found exactly. In addi¬ 
tion, we have developed systematic perturbation theory 
to account for corrections due to the short-range scat¬ 
tering. Concerning the transport properties, for weak 
interactions and magnetic field the magnetoresistance is 
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expected to be positive for a single Weyl node. (As com¬ 
pared to a negative magnetoresitance for current flow 
along the magnetic field when the realistic multivalley 
situation is considered.!^ However, based on the pre¬ 
sented analysis, one should not exclude a possibility of 
nonmonotonic magnetoresistance in systems with rela¬ 
tively strong interactions. 
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